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Direct photon spectra from central Au+Au- and Cu+Cu-collisions at ^/snn = 62.4, 130 and 
200 AGeV are calculated within the microscopic transport model UrQMD and a micro+macro 
hybrid model. In the latter approach, the high-density part of the transport evolution is replaced 
by an ideal 3+1-dimensional hydrodynamic calculation. We study the impact of viscosity and 
full local thermalization and compare the calculations to measurements obtained by the PHENIX 
collaboration. We find a reasonable agreement with the experimental data for calculations involving 
a Quark-Gluon-Plasma phase. 



I. INTRODUCTION 

Heavy-ion physics is widely used as a tool for the ex- 
ploration of the phase diagram of strongly interacting 
matter. In the collision of heavy nuclei, the nucleons 
may be compressed and heated sufficiently to create a 
new state of matter that consists of partonic degrees of 
freedom, the Quark-Gluon-Plasma (QGP) (T), Indeed, 
proposed signatures for the QGP, like strong jet quench- 
ing and large elliptic flow have been found by experiments 
at the Relativistic Heavy Ion Collider (BNL-RHIC) 

Inferring knowledge about the central regions of a 
heavy- ion collision is very difficult, since even if a plasma 
is created, its lifetime and size are beyond the experi- 
mental reach for direct observation, so we are limited to 
the study of particles that are emitted from the reaction 
zone. Unfortunately, first principle calculations of QCD- 
processes are only possible if all involved scales are much 
larger than the QCD-scale Aqcd ~ 0.2 GeV. However, 
in a heavy-ion collision, most particles have momenta 
comparable to Aqcd- Therefore, more phenomenologi- 
cal approaches are necessary to explore the bulk of the 
matter. 

While the abundance of hadronic particles that are 
produced in a heavy-ion collision are emitted at the end of 
the reaction and carry only indirect information from the 
early stages, electromagnetic probes allow for an undis- 
turbed view into all stages of the reaction. Photons and 
leptons escape the reaction zone without rescattering due 
to their very small cross-section, but for the same reason, 
their abundancies are rather low, compared to hadronic 
species 

Three different electromagnetic particle species arc 
currently being measured in heavy-ion experiments: 
single- and dielectrons, single- and dimuons and photons. 
Direct photons have the advantage that they are created 
in scatterings of the partonic or hadronic medium and 
are therefore directly coupled to the region of interac- 
tion. The leptons, however, are usually created in pairs, 
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either in the (initial state) Drell Yan process or by the de- 
cay of hadrons. In addition, one of the leptons might be 
a neutrino, which escapes observation. Since this process 
is governed by the weak interaction, the decay usually 
happens outside the fireball. Single leptons are therefore 
used to reconstruct weakly decaying heavy quarks, while 
the invariant mass distribution of dileptons can be used 
to extract spectral functions of vector mesons. 

Previous calculations of direct photons from transport 
theory include work with UrQMD by Dumitru et al. [Io| 
and Bauchle et al. and with HSD by Bratkovskaya 
et al. [l2T ] . Hydrodynamics has been used in many direct 
photon calculations, see e.g. p~3l - 420j . 

The extraction of the yield of photons from the fireball 
(direct photons) is hindered by a huge background of 
photons from hadronic decays outside the fireball, which 
is dominated by the ir°- and ^-decays. However, exper- 
imental techniques for the extraction of direct photon 
yields are well developed and allow to disentangle these 
late stage contributions from the scattering contribution. 
The experimental methods include a direct estimation of 
the background via invariant mass-analysis of the pho- 
tons [21], [22| , the analysis of interference patterns (using 
a Hanburry Brown- Twiss analysis) [23| and the extrapo- 
lation of the spectra of low-mass dileptons to the photon 
point [H. 

In this paper, we apply a previously established model 
for direct photon emission from hadronic and partonic 
sources and apply it to collision systems measured 
by the STAR and PHENIX collaborations at BNL-RHIC. 
In Section |nl we briefly introduce the model and the pa- 
rameters used for the present calculations, and in Sec- 
tion [HT]wc show the direct photon spectra obtained with 
our calculations as well as comparisons to the available 
data from the PHENIX collaboration Hi- 



ll. THE MODEL 

In the present work, direct photon spectra are calcu- 
lated in the framework of the microscopic Ultrarelativis- 
tic Quantum Molecular Dynamics (UrQMD) transport 
model (27l - [29| , using the hybrid option introduced in ver- 
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TABLE I: The critical energy densities for the mapping from 
hydrodynamics to transport theory for the various Equations 
of state, eo = 146 MeV/fm 3 is the nuclear ground state energy 
density. 

sion 3.3 US-HI [3|. While UrQMD itself is a hadronic 
transport model that includes only hadronic and string 
degrees of freedom and employs P YTHIA [33| for scatter- 
ings at high momentum transfer, the hybrid option allows 
to substitute the high-density part of the evolution by 
a 3+1-dimensional ideal hydrodynamic [3l[ description. 
In this part, other-than-hadronic degrees of freedom and 
phase transitions may be included. 

The inclusion of an intermediate phase into the model 
raises the need for two interfaces, to go from the particle- 
based description of the transport model to the density- 
based description of the hydrodynamic model and back 
again. 

The mapping from transport simulation to hydrody- 
namics is performed at i s tart = 0.6 fm. Here, the energy- 
density, baryon number-density and momentum densities 
are calculated from all particles at midrapidity. Particles 
with a rapidity \y\ > 2 are propagated in the cascade and 
do not interact with the bulk medium. 

The transition from hydrodynamics back to the cas- 
cade proceeds gradually, mapping the temperatures and 
chemical potentials to particles via the Cooper- Frye- 
formula [34[ when all cells in the same transverse slice 
(i.e. at the same position along the beam direction) have 
diluted below a critical energy density (see Table UJ. Af- 
ter the transition to the cascade, rescatterings and de- 
cays arc calculated in the well-known UrQMD model. 
For more detailed information on the hybrid model the 
reader is referred to [H, . 



A. Equations of State 

Three different Equations of State (EoS) are compared 
in this work. The effects of thcrmalization at the transi- 
tion from the initial stage cascade to hydrodynamics can 
be explored with the Hadron Gas- EoS (HG-EoS) [3(|. 
which has the same degrees of freedom as the transport 
phase. To investigate the effects of partonic matter and 
a phase transition, we use two different models for the 
EoS: The Chiral Equation of State x-EoS H3 has a 
cross-over phase transition to chirally restored and de- 
confined matter, while the Bag Model Equation of State 
BM-EoS [3l| has a first order phase transition to a Quark 
Gluon Plasma. In both EoS, the transition happens at 
around T c w 170 MeV. 



Centrality 


T s i op e [MeV] 


A [GeV~ 2 ] 


x7d.o.f. 


00%- 10% 


231. 9± 9.4 


2.39±0.67 


0.038 


00%-92% 


231. 4± 8.5 


0.41±0.11 


0.032 


10%-20% 


234.0±10.0 


1.26±0.37 


0.041 


20%-30% 


239.0±11.4 


0.56±0.18 


0.049 


30%-40% 


239.0±13.1 


0.27±0.10 


0.065 


40%-50% 


243.0±13.4 


0.12±0.04 


0.064 


50%-60% 


235.4± 8.8 


(5.64±1.43)-10~ 2 


0.032 


60%-92% 


250.5±11.8 


(6.91±2.08)-10" 3 


0.044 



TABLE II: Fit results for the low-p±-part (p± < 2.5 GeV) 
of the cascade calculations of Au+Au-collisions at ^/smn = 
200 GeV (see Fig. [TJ. The fit function is f(p±) = 



B. Photon emission sources 



Due to the small creation probability of direct photons, 
their emission is calculated perturbatively. I.e., the evo- 
lution of the underlying event remains unaltered by the 
emission of direct photons. 

The set of channels for direct photon production differ 
in the transport and hydrodynamic parts of the model. 
The most important channels, though, are common to 
both parts, namely 7T7t — > jp and irp — > jtt. Besides 
photon emission from the Quark-Gluon-Plasma, channels 
with strangeness are included in the hydrodynamic part. 
The corresponding rates for photon emission from each 
hydrodynamic cell are taken from Turbide et al. (l5j . In 
the transport part, additional processes including an 77- 
meson are included. The corresponding cross-sections 
have been calculated by Kapusta et al. [13| . 

Although Kapusta and Turbide use different La- 
grangians to derive their cross-sections and rates, ear- 
lier investigations (see [ll|) have shown that the thermal 
rates that can be extracted from Kapusta's cross-sections 
using this model agree very well with those parametrized 
by Turbide et al.. The same investigations have shown 
that the contributions of the hadronic processes that are 
not common to both models contribute about equally, 
but not significantly to the final spectra. The numerical 
implementation for direct photon emission is explained 
in detail in [ll| . 

At high transverse momenta, another source becomes 
important, namely the prompt contribution from hard 
scatterings of partons in the initial nuclei. The spec- 
tra predicted by NLO-pQCD calculations from Gor- 
don and Vogelsang J38| fit the experimental data from 
the PHENIX-collaboration [25[ rather well at high p±. 
Therefore, the pQCD contributions from (38|, scaled by 
the number of binary collisions (N co \\) are added to the 
soft photons calculated here. 
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FIG. 1: (Color Online) Comparison of the data from the 
PHENIX-collaboration [2f§| (black squares) to cascade cal- 
culations (red solid lines) for central to peripheral colli- 
sions. The green dash-dotted lines show the sum of pQCD- 
calculations [25], [38[ and the cascade contribution. For 
the most central collisions, 00-10% and 10-20%, the spec- 
tra from hybrid calculations with the BM-EoS plus pQCD- 
contribution are shown (violett dotted lines). 



III. RESULTS 

The comparison between direct photon spectra at low 
and intermediate transverse momentum p± from cas- 
cade calculations and data from the PHENIX collabo- 
ration (2f| for Au+Au-collisions at ^/.snn = 200 GeV is 
shown in Fig. [TJ One clearly observes that the hadronic 
transport model (full lines) docs not saturate the upper 
limits of the experimental data. In all centrality bins, the 
prompt photon yield is significantly larger than predicted 
by the hadronic cascade. The ratio between pQCD and 
hadronic contributions is fairly constant among the cen- 
trality bins. For comparison, Fig. [T] also shows the spec- 
tra obtained with the hybrid model using the Bag Model 
EoS (BM-EoS) for the two most central bins, 00-10% and 
10-20%, which agrees nicely with the data. Thermal fits 
to the low-pj_-parts of the cascade spectra show inverse 
slope parameters of T s | ope sa 235 MeV throughout the 
centrality bins, see Table HI1 

A more detailed exploration of the low-p_i_-part of the 
direct photon calculation is shown in Fig. [2j Here, the 
low-p^-data obtained by extrapolating the dilepton yield 
to zero invariant mass [26[ for central (00-20%) and mid- 
central (20-40%) is shown in comparison to cascade calcu- 
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FIG. 2: (Color Online) Comparison of data from the 
PHENIX-collaboration [26( (black squares) to cascade cacu- 
lations (red crosses) and hybrid-model calculations with HG- 
EoS (blue solid lines), x-EoS (orange dashed lines) and BM- 
EoS (violett dotted lines) for central (0-20 %) and mid-central 
(20-40 %) collisions. The contribution from initial pQCD- 
scatterings [26, [3£| have been added to all spectra. The spec- 
tra from central collisions have been scaled by a factor of 10 3 
to enhance readability. 



lations (red crosses) and hybrid calculations with hadron 
gas EoS (HG-EoS, solid blue lines), chiral EoS (%-EoS, 
dashed orange lines) and bag model EoS (BM-EoS, dot- 
ted violett lines) and prompt (pQCD) photon calcula- 
tions. All calculated spectra include the (A^i^-scalcd 
prompt photon contribution. 

In both centrality-bins, the direct photon spectra ob- 
tained with the BM-EoS and x _ EoS, which include a 
phase transition to a dcconfined state of matter, are sig- 
nificantly higher than the hadronic HG-EoS-calculations 
and agree with the measured data. 

A similar picture presents itself in Au+Au-collisions 
at lower incident energy ^/snn = 62.4 GeV and ^/snn = 
130 GeV, shown in the upper panels of Figure [3] The 
cascade calulations have been omitted from the Figure 
for clarity. 

Thermal fits to the spectra (see Table llllj) show in- 
verse slope parameters in the range from 233 < T s | ope < 
262 MeV, with the cascade calculations showing the 
smallest and the x _ EoS hybrid calculations showing the 
largest values of T 5 | ope . HG-EoS and BM-EoS calculations 
show similar inverse slope parameters. The integrated 
yield A is highest in BM-EoS hybrid calculations. The 
spectra from the hybrid calculations are rather similar 
for the different beam energies. 

Hybrid model calculations for central (0-20%) and mid- 
central (20-40%) Cu+Cu-collisions are shown in the lower 
panels of Figure [3] for all EoS. The thermal fits (see Ta- 
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P± [GeV] p ± [GoV] p ± [GoV] 

FIG. 3: (Color Online). Direct photon spectra calculated 
with the Hybrid model and HG-EoS (solid blue lines), \~ 
EoS (dashed orange lines) and BM-EoS (dotted violett lines) 
without prompt photon contribution. The left panels show 
calculations for ^/smn = 62.4 GeV, the middle panels show 
calculations for ^/snn = 130 GeV and the right panel shows 
calculations for ^/snn = 200 GeV. The upper panels show 
calculations for Au+Au-collisions, while the lower panel show 
calculations for Cu+Cu-collisions. In each panel, the upper 
curves are central collisions (00-20%) and the lower curves are 
mid-central collisions (20-40%). 

blc IIV|) again show no significant energy dependence of 
inverse slope parameter T 5 | p e or yield A. We observe a 
clear ordering of the total yield between the Equations 
of State, with yield from the BM-EoS calculations being 
higher than that of the x _ EoS, and both yields exceeding 
that of HG-EoS calculations. However, the inverse slope 
parameters are similar in HG-EoS and %-EoS calculations 
but significantly lower in BM-EoS calculations. 



IV. SUMMARY 

We examined the direct photon spectra obtained 
with a transport and a transport +hydrodynamics hy- 
brid model for collisions of Au+Au and Cu+Cu at en- 
ergies of ySi^ = 62.4, 130 and 200 GeV. We find that 
the hadronic models (transport model and hybrid model 
with Hadron Gas EoS) underpredict the data, while cal- 
culations with a dcconfined state of matter (hybrid model 





EoS 


centr. 


T slope [MeV] 


A [GeV -2 ] 


x 2 

d.o.f. 


200 


Transport 


0-20% 


232. 5± 9.8 


1.65±0.48 


0.041 


200 


HG-EoS 


0-20% 


246.7± 8.6 


3.63±0.83 


0.025 


200 


X-EoS 


0-20% 


261. 9± 8.7 


10.13±2.05 


0.020 


200 


BM-EoS 


0-20% 


251. 4± 9.7 


16.37±4.03 


0.029 


200 


Transport 


20-40% 


237.3±12.1 


0.38±0.13 


0.057 


200 


HG-EoS 


20-40% 


243.4± 8.3 


1.32±0.30 


0.025 


200 


X-EoS 


20-40% 


253.0± 8.0 


4.11±0.82 


0.020 


200 


BM-EoS 


20-40% 


240.6± 9.0 


7.61±1.90 


0.030 


130 


Transport 


0-20% 


232. 5± 9.1 


(9.87±2.67)* 


0.035 


130 


HG-EoS 


0-20% 


246.3± 8.5 


3.42±0.66 


0.024 


130 


X-EoS 


0-20% 


261. 2± 8.5 


9.67±1.93 


0.019 


130 


BM-EoS 


0-20% 


250.2± 9.6 


15.84±3.88 


0.039 


130 


Transport 


20-40% 


257.2±11.3 


(5.48±1.50) + 


0.036 


130 


HG-EoS 


20-40% 


242. 4± 7.6 


1.26±0.26 


0.021 


130 


X-EoS 


20-40% 


252.7± 7.9 


4.01±0.80 


0.019 


130 


BM-EoS 


20-40% 


240.6± 8.8 


7.46±1.82 


0.029 


62.4 


Transport 


0-20% 


242.1±13.5 


(5.29±1.95)* 


0.066 


62.4 


HG-EoS 


0-20% 


247.3± 8.1 


3.19±0.67 


0.022 


62.4 


X-EoS 


0-20% 


261. 8± 8.2 


9.24±1.78 


0.018 


62.4 


BM-EoS 


0-20% 


250.3± 9.5 


15.13±3.65 


0.028 


62.4 


Transport 


20-40% 


232. 8± 9.4 


(4.18±1.16) + 


0.038 


62.4 


HG-EoS 


20-40% 


245.8± 8.0 


1.21±0.26 


0.022 


62.4 


X-EoS 


20-40% 


253.9± 7.7 


3.82±0.73 


0.018 


62.4 


BM-EoS 


20-40% 


240.8± 8.6 


7.33±1.74 


0.028 






*: X10 


" 2 , +: X10" 


-3 





TABLE III: Fit results for the low-p ± -part (p± < 2.5 GeV) 
of the spectra from central (0-20%) and mid-central (20- 
40%) Au+Au-collisions. The fit function is f(p±) = 

A exp ^ — ~r~^ — ) ' ^ ne data are shown in Figure[2](for \J snn = 
200 GeV) and Figure |3] (^/snn = 62.4 GeV, upper left panel 
and yfstm = 130 GeV, upper central panel) . 



with Chiral or Bag Model EoS) fit the data much better. 

Thermal fits to the data show no significant beam en- 
ergy dependence on the spectra. The inverse slope pa- 
rameters obtained by fitting the low-transverse momen- 
tum part of the spectra are in the range of 227 < T s | 0pe < 
262 MeV, which is significantly above the expected tran- 
sition temperature to deconfined matter. 

Prompt photons from the initial early hard proton- 
proton scatterings are found to be a significant source 
of direct photon emission above p± = 3.5 GeV, if an 
EoS with phase transition is assumed, and is dominant 
throughout all p± if a purley hadronic scenario is as- 
sumed. 



V. OUTLOOK 

Future work with this model will include the extrac- 
tion of radial and elliptic flow parameters v\ and V2 for 
more differential analyses. Also, the influence of chang- 
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V snn 


EoS 


centr. 


T slope [MeV] 


A [GeV" 2 ] 


x 2 

d.o.f. 


200 


HG-EoS 


0-20% 


252. 0± 


9.6 


(4.84±1.38) x 


0.057 


200 


X-EoS 


0-20% 


251. 5± 


7.3 


1.77±0.39 


0.033 


200 


BM-EoS 


0-20% 


237. 7± 


7.8 


3.61±0.94 


0.047 


200 


HG-EoS 


20-40% 


254.6±13.2 


(1.61±0.62) x 


0.103 


200 


X-EoS 


20-40% 


242.9± 


7.0 


(7.25±1.63) x 


0.036 


200 


BM-EoS 


20-40% 


229. 2± 


7.5 


1.60±0.43 


0.051 


130 


HG-EoS 


0-20% 


250.0± 


9.3 


(4.78±1.35) x 


0.056 


130 


X-EoS 


0-20% 


250.9± 


7.1 


1.76±0.37 


0.031 


130 


BM-EoS 


0-20% 


238. 1± 


7.9 


3.56±0.93 


0.048 


130 


HG-EoS 


20-40% 


240. 4± 


7.7 


(1.99±0.50) x 


0.044 


130 


X-EoS 


20-40% 


242. 8± 


7.1 


(6.99±1.59) x 


0.036 


130 


BM-EoS 


20-40% 


228. 5± 


7.7 


1.58±0.44 


0.054 


62.4 


HG-EoS 


0-20% 


248. 2± 


7.7 


(4.71±l.ll) x 


0.039 


62.4 


X-EoS 


0-20% 


250. 2± 


6.8 


1.71±0.35 


0.029 


62.4 


BM-EoS 


0-20% 


236.8± 


7.4 


3.52±0.88 


0.044 


62.4 


HG-EoS 


20-40% 


242. 8± 


6.9 


(1.87±0.41) x 


0.034 


62.4 


X-EoS 


20-40% 


241. 7± 


6.3 


(6.71±1.37) x 


0.029 


62.4 


BM-EoS 


20-40% 


227.0± 


6.7 


1.62±0.40 


0.042 



x : Xl0 _ 



TABLE IV: Fit results for the low-p ± -part (p ± < 2.5 GeV) 
of the spectra from central (0-20%) and mid-central (20- 
40%) Cu+Cu-collisions. The fit function is f(p±) = 

A exp ( — Tffi- 1 — ) . The data are shown in Figure O lower pan- 

\ -'slope J 

els. 



ing the criteria for the transition between the transport- 
and hydrodynamic phases in the hybrid model will be 
examined in the future. 
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